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ABSTRACT 

Galaxy harassment has been proposed as a physical process that morphologically 
transforms low surface density disc galaxies into dwarf elliptical galaxies in clusters. It 
has been used to link the observed very different morphology of distant cluster galaxies 
(relatively more blue galaxies with 'disturbed' morphologies) with the relatively large 
numbers of dwarf elliptical galaxies found in nearby clusters. One prediction of the 
harassment model is that the remnant galaxies should lie on low surface brightness 
tidal streams or arcs. We demonstrate in this paper that we have an analysis method 
that is sensitive to the detection of arcs down to a surface brightness of 29 Bfj, and then 
use this method to search for arcs around 46 Virgo cluster dwarf elliptical galaxies. 
We find no evidence for tidal streams or arcs and consequently no evidence for galaxy 
harassment as a viable explanation for the relatively large numbers of dwarf galaxies 
found in the Virgo cluster. 

Key words: dust, extinction - dwarf galaxies - Virgo cluster. 



> 

in 
oo 

o 

l> 
o 
l> 
o 



1 INTRODUCTION 

Cold Dark Matter (CDM) hierarchical galaxy and large 
scale structure formation models predict many more small 
dark matter halos than can be associated with known dwarf 
galaxies (Moore et al., 1999a, Klypin et al., 1999). Methods 
of avoiding this discrepency within the models usually in- 
volve the global suppression of star formation in small halos, 
either by the prevention of infalling gas and/or the early ex- 
pulsion of gas (Dekel & Silk, 1986, Efstathiou, 1992). These 
two plausible physical mechanisms are critically challenged 
because the lack of large numbers of dwarf galaxies is en- 
vironmentally dependent. In the Local Group the ratio of 
dwarf (-14 < Mb < -10) to giant (-19 > Mb) galax- 
ies (the Dwarf to Giant Ratio or DGR) is about 5 (Mateo, 
1998). A similar value has been found by us for galaxies 
in other relatively low density environments (Roberts et al., 
2004a). The DGR is very different in the much higher density 
environment of the Virgo cluster, being about 20 (Sabatini 
et al., 2003, Roberts et al. 2004b)). There is good evidence to 
suggest that other dense environments also have relatively 
large numbers of dwarf galaxies (Kambas et al., 2000, Bern- 
stein et al., 1995). 

In a hierarchical model a cluster such as Virgo is con- 
structed from infalling galaxies, like the Milky Way and its 
companions, but the constituents of the cluster now are 
very different to that expectation. So, within the current 
paradigm, the cluster environment must have significantly 
affected its galaxy population in a way that does not affect 
galaxies in the field. This point has, of course, been known 



for many years because the centres of clusters are where we 
predominately find elliptical galaxies which, according to hi- 
erarchical models, are assembled from the merger of gas rich 
disc galaxies. Distant clusters also have a higher fraction of 
very blue galaxies compared to the field, again showing that 
the cluster environment is having an effect on the way galax- 
ies evolve (Butcher & Oemler, 1984, Dressier et al., 1994). 
What is still not clear though is just why there are so many 
more dwarf galaxies per giant in these rich environments. 
This has led to a number of ideas about how the cluster en- 
vironment may either make more dark matter halos become 
visible in the form of dwarf galaxies or how it may lead to 
the creation of dwarf galaxies. 

Galaxy 'squelching' is one proposed mechanism for 
making cluster dark matter halos more visible (TuUy et al., 
2002). The basic idea is that dwarf galaxies in a large over- 
density (that will eventually form a cluster like Virgo) form 
early, before re-ionisation, so that they can retain gas and 
form stars. Those outside of proto-clusters form later and the 
gas is then too hot to fall into small dark matter halos. TuUy 
et al. say that there is 'qualitative' agreement between their 
model and observation, but this was before the WMAP re- 
sult pushed the re-ionisation epoch back to z = 20 (Spergel 
et al., 2003), making the idea much less convincing. Babul 
& Rees (1992) have proposed that in the cluster environ- 
ment the expulsion of gas by dwarf galaxies is inhibited by 
the intra-cluster gas. This pressure confinement means that 
cluster dwarf galaxies are much more likely to retain gas 
to form subsequent populations of stars. Although a viable 
physical mechanism, we have shown (Sabatini et al., 2004) 



© 2003 RAS 



2 J. /. Davies et al. 



that this mechanism is only effective within the very central 
regions of a cluster like Virgo, where the intra-cluster gas 
density is high enough. But, this is also the region where 
dwarf galaxies would be totally destroyed by tidal forces 
(Sabatini ct al., 2004). Given also the competition between 
pressure confinement and ram pressure stripping it is not 
clear whether this is actually a viable mechanism at all even 
if dense gas were to exist outside of the cluster core. 

An alternative scenario is galaxy 'harassment'. This is 
a process that creates dwarf galaxies in the cluster environ- 
ment by the transformation of low surface density disc galax- 
ies. Numerical simulations by Moore et al., (1996, 1999b) 
indicate that the numerous high speed gravitational inter- 
actions that affect low density discs as they move through 
the cluster environment can remove the outer disc stars and 
'fatten up' the disc. The result is very much like the nu- 
merous dwarf elliptical (dE) galaxies that we see in clus- 
ters. Moore et al., have proposed that the 'excess' dE galaxy 
population of clusters like Virgo is duo to this morphologi- 
cal transformation of low surface density disc galaxies. It is 
clear that tidal interactions of some form or another must 
play an important role in cluster galaxy evolution. Observa- 
tions of distant clusters clearly show many apparently tidally 
disrupted galaxies (Dressier et al., 1994). The origin of intra- 
cluster stars (Ferguson et al., 1998) and planetary nebulae 
(Feldmeier et al., 1998) are presumably the result of such 
interactions. In summary there is good evidence for tidal in- 
teractions (Moore ct al. 1998), but is the evidence for mor- 
phological transformation so strong ? 

In this paper we look at one of the predictions of the 
harassment model to scrutinise it in the light of now observa- 
tional data and data reduction techniques. The harassment 
(morphogical transformation) model predicts that stars re- 
moved from the harassed galaxies will orbit in 'prominent' 
narrow streams that follow the orbital path of the galaxy. 
The galaxy itself is a brightening on this stellar stream. A 
rough calculation predicts that an annulus of stars tidally 
stripped from a Mb = — 20 (Moore et al., 1996) galaxy dis- 
tributed around a radius of 1 Mpc (of order the size of a 
galaxy cluster) and thickness 2 kpc (of order the size of a 
dwarf galaxy) would have a surface brightness of about 27.5 
Bjj, (Bjjt - blue magnitudes per sq arc sec). Features like this, 
that almost certainly have a tidal origin, have been found in 
a number of clusters. For example Calcaneo-Roldan et al., 
(2000) describe CCD observations of a low surface bright- 
ness arc in the Centaurus cluster (originally discovered by 
David Malin). The arc has a length of about 100 kpc and a 
width of 2 kpc, its surface brightness is 27.8 S//. Trentham 
& Mobaslicr (1998) describe an arc in the Coma cluster that 
is about 80 kpc long and has a surface brightness of order 
26.5 Bn- Gregg and West (1998) have identified other possi- 
ble tidal features in the Coma cluster core that have surface 
brightnesses of order 27.5 B/x. Though convincingly tidal in 
origin none of them are convincingly assiciatcd with a dE 
galaxy. But, the observations do indicate that if we want to 
search for arcs as signatures of harassment we need to be 
sensitive to features that have surface brightnesses as low as 
28 B/x. 

Should wo find such features in a cluster like Virgo ? 
According to Calcaneo-Roldan et al., (2000) these tidal fea- 
tures should persist for a few cluster crossing times, which 
for Virgo is a few billion years (Trentham et al., 2001). 



Moore et al., (1996) say that the dramatic change in clus- 
ter galaxy morphology due to galaxy harassment has oc- 
cured from z « 0.4 to today. This is a period of 4.6 billion 
years (standard flat Qm = 0.3, Ao = 0.7 and Ho = 70 km 
s^"*^ Mpc^^ cosmology). Thus our test, as described below, 
should be sensitive to galaxies that have been harassed over 
the latter half of the morphological transformation epoch. 



2 DATA AND DATA PROCESSING 

We have used B band CCD data from the Issac Newton 
Telescope (INT) Wide Field Camera survey of the Virgo 
cluster to search for evidence of tidal streams. These data 
arc fully described in Sabatini et al., (2003) and Roberts et 
al., (2004a). As the basis of our search for tidal features we 
have used 38 galaxies listed as dE in Roberts et al., (2004b). 
Of these galaxies 24 (63%) are also listed in the Virgo Clus- 
ter Catalogue (VCC) (Binggeh et al. 1985) as cluster dE 
galaxies. The typical one sigma sky noise for these frames is 
26 Bjj, some 2 magnitudes brighter than required to detect 
tidal arcs. 

To increase the surface brightness limit to 28 B/i we 
have smoothed the data. Smoothing can be problematic 
when looking for Low Surface Brightness (LSB) features as 
excessive smoothing, for example with a Gaussian or tophat 
function, makes all objects appear to be LSB. Ideally, one 
would like to smooth the 'sky' (those regions that contain no 
resolved foreground or background objects) without smear- 
ing out light from the resolved objects. One way of doing this 
is to apply the Haar transform. In its simplest form it can 
be used to smooth over a 2 x 2 pixel area. Instead of storing 
in memory the 4 pixel values 4 different numbers arc stored, 
from which the original image can be reconstructed - the 
mean value of the four pixels, a measure of the x-gradient, 
a measure of the y gradient and a measure of the curvature 
(for an example see Aboufadel and Sclicker, 1999). If the 
measures of the gradients and curvatures are as expected 
from a separate measurement of the background pixel-to- 
pixel fluctuations (they are just sums of differences of pixel 
values) then they can be set to zero or reduced. Where the 
gradients are high (when the pixels are part of a resolved 
object) the differences can be retained. After transforming 
back to pbcel values, the sky has been smoothed without 
spreading out the light from resolved objects. To generalise, 
we form the orthonormal matrix T which by definition has 
unit magnitude and satisfies = T^^. If the 2" x2" matrix 
X contains the pixel values then we can form the matrix 

Y = TXT"^ (1) 

in the 2x2 case the matrix Y contains the mean value 
and gradients described above. The matrix Y can now be 
filtered as described above and then transformed back using 

X = T'^YT (2) 

For the purpose of this experiment we want to filter on 
kpc scales (width of the tidal features) and so we bin the 
data into 1 arc sec pixels and then filter with a matrix T of 
dimension 16 x 16 (16 arc sec is ~1 kpc at the distance of 
Virgo, 16 Mpc). In the case of a 16 x 16 matrix, Y contains 
one element that is the sum of the pixel values and the 
others are the sums of the differences between pixels - so 
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Figure 1. Data with simulated arc (left) and after Haar smoothing (right), a) No arcs added, b) Diagonal line (top left to bottom right), 
c) Diagonal line (bottom left to top right), d)Horizontal line, e) Vertical line. Each is displayed between la below sky level to 3cr above 
sky level (cr is the sky standard deviation). 



again we can see if this is consistent with being sky noise 
and filter as appropriate. An example of the effect this has on 
a data frame is shown in Fig. la, it is clear that the brighter 
compact objects have not been smoothed. The filtering does 
tend to break the image up into 16 x 16 blocks, but we 



will show below that it still enables us to detect features at 
surface brightness levels of 28 Bfi. 

To test the method we have added 'arcs' of surface 
brightness (with random noise) of 28 Bjj, to data frames 
from the INT. We have extracted regions of 256 x 256 1 
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arc sec pixels centred on our sample of dE galaxies. This 
corresponds to about 20 kpc at the distance of the Virgo 
cluster (same linear size as the arc described in Conselice 
and Gallagher, 1999). We have chosen this distance because 
we want the 'arcs' to be almost linear over the region of the 
data we use. An 'arc' in the plane of the sky of radius 1 Mpc 
would only deviate from a straight line by one pixel over this 
distance. The arc described in Calcaneo-Roldan et al., 2000, 
deviates by 3-4 arc sec over ICQ arc sec - so we are expecting 
almost linear features. Our method (see below) would only 
then not work if we were unlucky enough to be viewing the 
galaxy close to the tangent point (orbital velocity projected 
along the line of sight) of an almost edge-on arc (orbit in- 
clined at ~ 90° to the plane of the sky) . The arc would then 
more resemble a 'tail' as it would extend to one side only 
of the galaxy. This will occur when the galaxy is less than 
about one data frame size away from the tangential point. 
For a 1 Mpc radius orbit this will amount to about 10% of 
the time. So, less than 10% of our target arcs might go un- 
detected because of the detection method we are using (see 
below). This is an upper limit because not all orbits will 
be sufficiently inclined. We also expect that the 'arcs' will 
both be brighter and stable over longer periods closer to the 
galaxy. In Fig. l(b-e). we show on the left a random piece 
of sky with faint 28 B/i 2 kpc wide arcs added. On the right 
we show the result of Haar smoothing. 



Looking at Fig. 1. it is not really clear that the Haar 
smoothing has helped very much in making the tidal features 
more visible. Below we will show that it has, we just need 
to quantify the effect. Given that our frames are centred 
on a dE galaxy and that we expect the arcs to be straight 
lines in our images we can now quantify our detections. After 
numerous trials, we consider only those pixels (in the filtered 
image) that have intensity values that lie between the sky 
background value and the sky background plus Scr (a is the 
background standard deviation in the filtered frame). We 
then calculate the angle that a vector connecting the pixel 
to the centre of the frame (position of the galaxy) makes 
with the horizontal. In Fig. 2. we show the distribution of 
angles derived from the simulated data shown in Fig. 1. Fig. 
2a) is for the smoothed data frame with no 'arc' added (a 
random area of sky) - no significant peak (see below) is seen. 
Fig 2b) is for a frame that has an 'arc' added, but has not 
been smoothed - again no significant peak is seen. Figs 2c-f. 
each shows a significant peak at the angle of the 'arcs' shown 
in fig. 1. We have quantified the significance of the detection 
by calculating the standard deviation of the distribution of 
angles (excluding the peak value) and comparing this to the 
number in the peak minus the mean value (the Signal-to- 
Noise (57V)). For the histograms shown in Fig 2a and b we 
find SN < 2, for those in Fig 2c-f SN > 3. To summarise, 
we make no detections in a frame on a random piece of sky 
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(checked with 20 other random pieces of sky) or on frames 
that have had 'arcs' added, but have not been smoothed. All 
simulated arcs are detected with a SN > 3 when smoothed. 
We conclude that on real data frames a value of SN < 2 
corresponds to a non-dctcction while SN > 3 is a definite 
detection. For 2 < SN < 3 the question remains open. 

We have applied the above detection method to our 
sample of 39 dE galaxies from the Virgo cluster. The result 
is that 26 (67%) had SN < 2 showing no sign of any tidal 
features. There were no galaxies with SN > 3, and therefore 
no definite detections, although 13 (33%) had 2 < SN < 3 
and so still remain possible detections. 



3 MOVING TO LOWER SURFACE 
BRIGHTNESSES 

One obvious critisism of the above technique is that the 
simulated axes were added to data that was already reduced. 
This was done because the data is pipeline processed by the 

Cambridge Astronomical Survey Unit and then the reduced 
data is accessed from their archive. To go both deeper in 
surface brightness and to simulate the tidal streams in a 
more realistic way we have obtained additional observations 
of 7 galaxies from the Virgo cluster sample (Sabatini et al., 
2003). 

In April 2004 we used the DOLORES instrument on the 
Galileo Telescope, La Palma (TNG) to obtain 50 minute 
exposures for each of our targets. The data was reduced 
in the normal way using the IRAF package and calibrated 
from photometric standards observed at regular intervals 
throughout the night (photometric accuracy to ±0.1 mag- 
nitudes). We chose to use the R band because with this 
instrumental set up we could go deeper in a given exposure 
time at R than at B. Typical early type galaxy colours are 
1.0 < {B — R) < 1.5 which we will use to compare the B 
band observations of the previous section with the R band 
observations discussed here. 

To test the procedure described in the previous section 
we have added simulated arcs to the raw data (each image 
is made up of 5 separate 10 min exposures) and then pro- 
cessed the frame as before. Below we will show that an arc of 
surface brightness 27.5 Rp (28.5-29 B/x) is easily detectable. 
This is 0.5 to 1.0 magnitudes deeper than the observations 
described in the previous section. In Fig. 3 we show the data 
frame (before smoothing) and the distribution of angles that 
low surface brightness pixels make with the origin. Looking 
carefully, the simulated stream can just be made out as a di- 
agonal line across the CCD from bottom left to top right. In 
the histogram an arc at an angle of 45 degrees is clearly seen 
at SN = 4.4. It is quite clear that arcs at this surface bright- 
ness level are not removed by the data processing. Given the 
strong signal obtained for an arc at a surface brightness of 
27.5 R/i we are confident that our procedure is capable of 
detecting arcs in this data set to an equivalent level of 29 
B/x. 

We have applied the above procedure to the 7 galaxies in 
our sample. None of the galaxies have a detection at SN > 3; 
one has 2 < SN < 3. There is no evidence for any low surface 
brightness arcs. 



4 DISCUSSION 

Although the original numerical simulations of galaxy ha- 
rassment were carried out with a cluster designed to be like 
Coma, rather than Virgo, the conclusion was that this was 
a viable method for creating the excess of dwarf galaxies 
found in clusters in general (Moore et al., 1999). The clear 
result from this paper is that we find little or no evidence 
for any extended tidal streams, with surface brightnesses less 
than 29 B/u, associated with our target galaxies that would 
be the signature of morphological transformation by galaxy 
harassment. There are a number of possible explanations. 
The most obvious is that they do not exist (discussed in 
more detail below). The second is that they are too faint. 
This could arise because we have over estimated the stellar 
luminosity of the progenitor galaxies (the arcs are at much 
lower surface brightnesses) or because the tidal material was 
pulled out of the galaxies a long time ago (more than a few 
billion years) and has since dispersed. We cannot rule out 
either of these two, but we believe that there are other rea- 
sons for being skeptical about the harassment model - the 
morphological transformation of LSB disc galaxies into dE 
galaxies: 

(i) There needs to be a large population of low surface 

brightness disc galaxies that reside in the Universe at red- 
shifts of order z = 0.4, awaiting morphological transforma- 
tion. 

(ii) Moore et al. (1996) quote a number of articles that 
clearly demonstrate the Butcher-Oemler effect - that is that 

the fraction of galaxies with colours similar to late type 
galaxies is higher in distant clusters compared to nearby 
clusters. They then associate this with the morphological 
transformation of galaxies. It is clear that galaxies are un- 
dergoing enhanced star formation and can be morphologi- 
cally disturbed in the cluster environment, but the Butcher- 
Oemler effect is not evidence for morphological transforma- 
tion. 

(iii) There are many dE galaxies that are just too small 
to have been re-shaped by harassment. If the cluster dwarf 
galaxy population originated from a tidally truncated pop- 
ulation of larger field galaxies then their tidal radii should 
be « ^ fi^££i/ -vvhere Rc is the cluster core radius (« 0.5 

Mpc for Virgo), adwarf (~ 10 km s~^) is the internal veloc- 
ity dispersion of the dwarf and aciust (~ 700 km s~^) is the 
velocity dispersion of the cluster. The smallest tidal radius 
is of order 7 kpc. This is much larger than most of the dE 
galaxies we have in our Virgo cluster sample (radii of order 
10 arc sec, less than 1 kpc). 

We have previously proposed that dwarf galaxies are 
rather robust objects in the cluster environment because 

they have large mass-to-light ratios (Sabatini et al., 2004). 
However, in clusters like Virgo they do undergo accelerated 
evolution (star formation) due to tidal interactions in the 
cluster (see for example Henriksen & Byrd, 1996) which 
eventually exhausts their gas supply. Thus the cluster en- 
vironment is one that reveals its dark matter lialos to us. 
So, are cluster dwarf galaxies the embers of a once much 
brighter fire or are they ignited for the first time in the clus- 
ter environment ? 
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Figure 2. Histograms of the angles that significant pixels make with the frame centre pixel, a) The smoothed data frame with no arcs 
added, b) Unsmoothed frame with an arc added, c) Smoothed with arc at an angle of -45 deg, d) Smoothed with an arc at an angle of 
45 deg, e) Smoothed with arc at angle of deg, f) smoothed with arc at an angle of 90 deg. 
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show the histogram of the angles that significant pixels make with the frame centre. There is a strong signal at 45 deg. 
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